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^ ■ Abstract 

■ This article investigates asymmetric brane-world scenarios in the limit when the bulk 

(«~^ I gravity is negligible. We show that, even when true self gravity is negligible, local mass 

\^ ■ concentrations will be subject to a mutual attraction force which simulates the effect of 

I Newtonian gravity in the non-relativistic limit. Cosmological and also post-Newtonian con- 

straints are examined. 

:^ ■ 
"o ■ 

O^' 1 Introduction 

u 

^P' Gravity has long been considered to be one of four fundamental forces. Yet it is very different in 

its nature from the others and, although particle theories have been devised to unify the other 
^ , three, it has proved surprisingly difficult to unify all four. One possible way to cut this Gordian 

^ I knot is to consider gravity not as a fundamental force but as a derived one In this paper 

we consider one possible way to realise this within the context of higher dimensional theories 
by explaining gravity as an apparent effect due to the acceleration of a brane-world universe 
through a five dimensional spacetime. 

Following an original suggestion by Randall and Sundrum many authors have considered 
5-dimensional scenarios in which effectively 4-dimensional (inverse square law) gravity is induced 
on a 3-brane hypersurface representing our observed 4-dimensional spacetime due to the effective 
confinement of the 5-dimensional gravitational field by a strong extrinsic curvature. The kind 
of cosmological scenario that was originally proposed was characterised by a postulate of Z2 
reflection symmetry about the 3-brane hypersurface, in which case the effective 4-dimensional 
Newton constant Gj^j was found to be proportional to the square of the analogous 5-dimensional 
gravitational coupling constant Gj^j with a coefficient proportional to the asymptotic limit value 
%^ of the brane tension. A more general relation 
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has been obtained more recently IH allowing for the possible presence of a force density term /, 
such as might arise naturally from a gauge four form coupling [^. 

Moreover, for scenarios of this more general reflection symmetry violating kind |5|, ^, |7| it has 
subsequently been shown ^ that it is still possible to recover an effective 4-dimensional Einstein 
equation for local perturbations on the brane under the same kind of conditions as in the earlier 



demonstration for the more specialised reflection symmetric case by Shiromizu et al [10|. The 
applicability of the effective 4-dimensional Einstein equations with the coupling constant 

given by (|l|) is restricted to the weak source limit (relevant for the present epoch but not for 
earlier stages of cosmological evolution during which gravity would have behaved in a "non- 
conventional" manner Q) and to the case where the background Weyl tensor contributions, 
representing the effect of incoming gravitational waves or black holes in the five dimensional 
"bulk" , are sufficiently small. 

The purpose of the present work is to consider a very different kind of limit, involving a 
brane- world scenario in which an inverse square law for gravity is recovered, not by an effective 
confinement mechanism due to strong extrinsic curvature, but simply by having a sufficiently 
small five dimensional gravitational coupling constant for genuine gravitation to be entirely neg- 
ligible. In this weak gravity limit our previous formula (||) will cease to be applicable otherwise 
an, unacceptable, negative value for Gj^j would be obtained: in the limit considered here, the 
derivation of (||) for localised perturbations |^] will break down because the Weyl contribution 
(on whose negligibility its validity depended) will become indeterminate due to Gj^j terms in the 
denominator. In order to deal with this limit, which is singular from the previous point of view, 
we need to start again with a fresh approach. 

The preceding work Q shows that the worldsheet geometry, and hence the apparent gravity 
on the brane, has three origins. There is a first part coming from the geometry of the "bulk" 
background, a second part arising from the discontinuity of the extrinsic curvature form, X^j^, 
across the worldsheet, and a third part coming from its mean value. In the reflection symmetric 
case that is most commonly considered, the latter effect is absent (that is, there is a vanishing 
mean {K^^) = 0). The aim of the present article is to investigate the opposite extreme limiting 
case characterised by the absence of gravitational perturbations and, in particular, of the grav- 
itationally engendered discontinuity of the extrinsic curvature (that is, we will have [K^y\ = 0) 
so that it is the mean curvature, identifiable in this case with the value on either side, which 
provides the dominant contribution. 

What will be shown here is that, when five dimensional gravity is absent or negligible, an 
artificial, and automatically positive, 4-dimensional (inverse square law) gravitation effect will 
still be present whenever reflection symmetry is broken, with an effective coupling "constant" 
that will be proportional to the rest frame time component, K^^ say, of the extrinsic curvature 
of the averaged background, according to a formula which, instead of (|l|) will have the form) 

However although it mimics Newtonian gravitation in the non-relativistic limit, the artiflcial 
gravitation occurring in the simple type of 5-dimensional 3-brane scenario considered here will 
produce post-Newtonian effects that will, in general, differ from those of the standard Einstein 
theory, typically in a parameter range beyond what is observationally admissible. 

Although it will be shown that it can be contrived in such a way as to satisfy the cosmological 
nucleosynthesis constraint and instantaneously mimicking the standard Einstein theory in the 
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immediate spacetime neighbourhood of our solar system, the necessary fine tuning would seem 
to be very artificial. Furthermore, a gravitational "constant" that is proportional to the, possibly 
variable, extrinsic curvature of the average background can hardly be compatible with the very 
severe observational limits on deviations from constancy of gravitational coupling. 

These considerations lead us to conclude that in its present form the simplest model con- 
sidered here does not by itself provide a viable alternative theory of gravity. However — quite 
apart from the possibility that a viable alternative may be provided by some generalisation — 
the simulated worldsheet gravity phenomenon considered here will still be of physical interest 
and potential practical relevance as a mechanism that can provide small but not necessarily 
negligible additions to effects of genuinely gravitational origin (in much the same way as, from 
a laboratory point of view, the genuine gravitational attraction of our Earth is modified by the 
centrifugal effect of its rotation). 

The plan of this article is as follows. The necessary worldsheet perturbation formalism for 
the general case of a codimension one brane of arbitrary worldsheet dimension p is developed 
in section |2[ The relevant non-gravitating kind of minimally coupled Q homogeneous isotropic 
cosmological background to which it will be applied is presented next in section |3|, where the 
first integrated (Friedmann type) form of the dynamical equations is obtained in terms of an 
adjustable constant of integration (whose interpretation as the energy with respect to the sta- 
tionary geometry of the external "bulk" is described in the appendix) . The heart of this article 
is in section Q where it is shown how local perturbations will generically be subject to an at- 
traction of Newtonian gravitational in the non-relativistic limit, and where the more restrictive 
conditions for agreement with the linearised Einstein equations are derived. The remainder 
of the discussion is restricted to the cosmologically interesting case of a hypermembrane with 
worldsheet dimension p = A, starting in section ^ where the conditions for compatibility with 
the standard description of nucleosynthesis are derived. In conclusion, the difficulty of reconcil- 
ing these global conditions with what is required for matching the parameters governing local 
gravitational effects is discussed in section ^ 

2 Lagrangian worldsheet perturbations 

As a technical prerequisite for the perturbation analysis that will be involved, we start by 
considering the effect of the brane worldsheet perturbation, which we suppose to be generated 
by an infinitesimal displacement vector field say, on the first and second fundamental forms 
of the worldsheet, which correspond to tensors denoted by 7^^ and iC^jy using the notation of 
our previous work In terms of the background spacetime metric g^,^ and the unit normal 
(^^ = V^Ci where C measures orthogonal distance from the brane (on the side chosen to be 
considered as positive) the first fundamental tensor is defined simply by 




(3) 



(4) 



the tensor corresponding to the second fundamental form is given[] by 



(5) 



^Wc arc using the MTW sign conventions 
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and will automatically have the Weingarten symmetry property K^y = Ky^ as a worldsheet 
integrability condition. This tensor is also obtainable from the general second fundamental 
tensor Kfj_yP that was used in an earlier brane perturbation analysis by K^^ = K^^f^p. 

Since our purpose here is to consider the limit in which genuine gravitation is negligible, the 
geometry will be globally unaffected by the perturbation, so the Eulerian metric perturbation 
hfiu = ^^dfiu that was allowed for in the earlier analysis can simply be set to zero. Thus 
locally, with respect to coordinates that are comoving with respect to the displacement (so that 
the brane worldsheet retains the same coordinate locus, C = 0) the only change in the metric 
tensor will be the purely Lagrangian variation given by the familiar Lie differentiation formula 

Sl^9^iu = 2V(^C^) , (6) 

using round brackets to indicate index symmetrisation. Preservation of the locus of the world- 
sheet implies preservation of the direction of the covariant normal C^, but the change @ in 
the metric will entail a corresponding change in amplitude to preserve the unit normalisation 
property C^lC^ = 1; so it can be seen that the ensuing Lagrangian variation of the unit normal 
will be given by 

s^Cf. = Wp'^''e- (7) 

This formula can be used to obtain a direct verification of the corresponding Lagrangian variation 
formula for the first fundamental tensor, namely 

6,^" = -^iJ^V'^k" , (8) 

which is actually valid, not just for the hypermembrane (codimension 1) case under consideration 
here, but quite generally for a hyperstring (the case with codimension 2, as in a 6-dimensional 



world brane scenario |14]) and cases with even higher codimension. 



Using the explicit expression 

"^f^C = C:, + r/, C , r/^ = g'''' {q.m - i^mp,.) , (9) 

for the covariant derivative that is involved in (|5|) and the well known fact that the Lie derivative 

^lVp = g""" (v(/,ff^), - \vj,g^.^ , (10) 
of the Christoffel connection can be expressed in terms of the background Riemann tensor TZfj,upa 



as 



Together with (0) and (^), this enables us to calculate the corresponding variation of the extrinsic 



curvature as 



d^Kf"" = Cp¥''kY\ [v^'^v^^e + n'^p^^c) + 4.Kyv''^c ■ (12) 



In the language of |1S], this can be obtained from the formula Kp^ = Kp^yfQp using (0) and 
the general expression for 6^Kp_yP . If one makes the convenient choice of an orthogonal gauge, 
setting 

e = s,e, (13) 
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it can be seen that ( |12|) will be expressible directly in terms of the scalar displacement amplitude 
^ as 



In what follows we shall be particularly interested in the trace K = KJ^, for which the cor- 
responding variation can be seen to be expressible in terms of the background Ricci tensor, 



5^K = V^'V^i + {K^'''K^,+ ll'^n^,)^. (15) 



This formula leads to a noteworthy cancellation when used in conjunction with the contrac- 
tion of the well known Gauss identity (see, for example, the preceding analysis |jl^) according 
to which the Ricci scalar R of the induced metric on the hypersurface will be given by 

R = K^- Rf'^K^^ + l^^^'n^u - n . (16) 

The preceding formulae are valid regardless of the dimension, p say, of the worldsheet under 
consideration, the only restriction being that it should have codimension one. However we must 
take account of this dimension in the next step, which depends on separating the trace-free part 
of the extrinsic curvature tensor, 

C^, = K^,-^K^^^, C = 0. (17) 

and also the trace-free part of the background Ricci tensor, 

S^,, = n^u - -^T^Qi^v , 5; = , (18) 

which vanishes whenever the background satisfies Einstein type vacuum equations, with or 
without a cosmological constant. What is remarkable is that the form of (ITsI ) is such that the 
identity (^) can be used to eliminate not just the extrinsic curvature trace K, but also the 
background Ricci trace scalar, IZ, in favour just of the intrinsic Ricci trace scalar i?, and hence 

5^K = WV^ i + ^ {r + p C^^'C^, -{p + l)Y"S^u) . (19) 

This generalises the results by Garriga and Vilenkin [^]. The particular case p = 2 of this 
identity applies to any ordinary 2-surface in 3-dimensional space, while the case p = 2> applies to 
the worldsheet of an ordinary membrane, including the case of a cosmological domain boundary, 
in 4-dimensional spacetime. What we shall be concerned with here is the case p = 4 of a 3-brane 
in a 5-dimensional spacetime background. In this latter context, the significant feature of (^) 
is that in typical brane-world scenarios the quantities in the bracket on the right are all zero or 
tend to zero at late cosmological times, unlike the cancelled traces K and 7^, whose amplitudes 
remain very large. 

3 Homogeneous cosmological reference configuration 

The previous analysis will allow us to study the Newtonian limit of the theory. First, we shall 
study homogeneous and isotropic reference configuration which will also be of interest for the 
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cosmological solutions. Our purpose now is thus to consider the case of a hyper membrane 
governed by an equation of motion taking the standard form 

T^'^i^/..=7, (20) 

that apphes ||l^ in the hmit when self gravitation is negligible so that the background geometry 
can be considered to be smooth. In this dynamical equation, / is the orthogonal component of 
the external force density, and T^^ is the surface stress energy density tensor, which will satisfy 
an intrinsic pseudo-conservation law of the form 

i.p%r'"' = f), (21) 

provided the external force density has no tangential component. We shall restrict our attention 
to the kind of force produced by minimal gauge coupling for which this orthogonality 

condition will automatically be satisfied, and for which the force density / will simply be a 
constant. 

We wish to study perturbations of a homogeneous isotropic reference configuration in a 
background satisfying Einstein type vacuum equations of the standard form 

5^. = . (22) 

As discussed in the preceding work , application of the relevant generalisation of the well 
known 4-dimensional Birkhoff theorem ||l^ allows us to conclude that the background metric 
will have the static form 

^s^ = r^Ae + '^-Vdt\ (23) 

where di"^ is the positive-definite space metric of a (p — l)-dimensional sphere, plane, or anti- 
sphere with constant curvature, k say, respectively positive, zero, or negative, and V is a function 
only of r. The location of the brane worldsheet will be given by an expression of the form r = a 
where a is a function only of t. There will be a equivalent brane based formulation 

ds2 = r2 + dC^ - i/^dr^ (24) 

in which the worldsheet locus is given simply by C = 0, and in which r represents proper time 
on this locus, so that v = \ there. Using a dot to denote partial differentiation with respect to 
r, and a dash for partial differentiation with respect to C) the rate of change of the scale factor 
a will satisfy the conditions 

= r'2 - V and = r'^/v'^ . (25) 

In such a system, the first fundamental form ^ specifying the induced metric on the worldsheet, 
is 

ds^ = dx^' dx" = dt^ - dr^ , (26) 
and the corresponding expression for second fundamental form (^) is 

K^^ dx^ dx'' = -ar' di^ + i/' dr^ , (27) 

which can also be expressed as 

Kf,^ = -—7^^ + f^' ~ ^M^'^ ■ (28) 
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in terms of the preferred timelike unit vector = "S/yT that specifies the cosmological reference 
frame. In terms of this unit vector, the stress-energy momentum tensor will have the isotropic 
perfect fluid form 

T^"" = -TY"" + {U- T)u^'if , (29) 

where lA is the total rest frame energy density and T is the total brane tension. It can be 
deduced from (^) that the dynamical equation (pO|) takes the form 

Uv' + {p- 1)T- = 7 (30) 
a 

from which we deduce that the extrinsic curvature scalar is given by 

while the trace-free part will have the form 

C^^v = C Qt^i/ + Uf,u^ , (32) 
so that the scalar in ( p^ is explicitly given by 

Ct'T...,. = iu' - - 

p — 1 

Due to homogeneity, the non-trivial content of the pseudo-conservation law (^l|) reduces to 
the simple energy equation 

U+ip-l)iU-T)- = 0. (34) 
a 

Using (^) and the relation r' = v'a obtained from (|2^), the dynamical equation (|30|) can be 
rewritten as 



C= = ^C-a^^ = ( 1/' - - I . (33) 



Since, as remarked above, the minimal coupling postulate automatically ensures ^ that / is 
constant, this equation will be immediately integrable to give 

where is a constant of integration, which can be interpreted as being proportional to the 
global energy of the system. Since the concept of energy normally involves reference to some 
stationary background reference system so as to be well defined, it can not be used in the usual 
global cosmological context of an expanding universe model. It will, however, be shown in the 
appendix how a bulk background which is static — due to the generalised Birkhoff theorem |P 
— will allow us to define this energy in an unambiguous manner for minimally coupled scenarios 
of the kind considered here. Using (30) again, we obtain the corresponding formula 
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which completes the evaluation of the coefficients in the expression for the second funda- 
mental form. Therefore, the extrinsic curvature scalar is given by 



and the corresponding expression for the scalar coefficient C in (32) and (^) is 

As in the preceding work ||3|, ^ on the self gravitating case, we adopt the usual supposition 
that the hypersurface stress energy tensor T^'^ in (^) has the form 

Yf^u ^j-^-t^u ^ (40) 

in which t^'^ is the part attributable to matter fields on the brane and ^ is a constant inter- 
pretable as the "bare" brane tension meaning the limit to which the tension and energy density 
tend in the infinitely distended limit. In the homogeneous case characterised by (^) the material 
contribution will be given in terms of average cosmological density and pressure functions 

-p = U-T^, P = T^-T, (41) 

by 

t"" = P 7^^^ + (7, + P) n^u'^ . (42) 

In the large expansion limit a — > oo that is relevant at late cosmological times the material 
contributions will tend to zero, p ^ and P — s- 0, so that lA ^ and T ^ which implies 
that we shall have _ 

K^K^^-^, C^O, (43) 

oo 

which will be discussed in more details in the following sections. 

4 Perturbed Configuration 

We now turn our attention to an "actual" configuration in which the brane worldsheet deviates 
from that of the homogeneous reference configuration described above by a small displacement 
with amplitude ^ of the kind described in section ^. Using a tilde to distinguish perturbed 
quantities from their analogues in the reference configuration, we take the first and second 
fundamental tensors of the perturbed worldsheet to be given, using the notation of ( |Tl] ) and 
(0), by 

l^lu = li,u + Kltiu > (44) 

and 

K^, = K^, + 5^K^,. (45) 
The corresponding perturbed version of the equation of motion (20) will have the form 

f'"'Kf,u=7, (46) 
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with a perturbed surface energy-momentum tensor of the form 



T^"' = -T^Y"" + , (47) 

in which the constants 1^ and / are just the same as in the reference configuration. 

By subtracting (^) from its unperturbed analogue ( pO|) we are left with a dynamical source 
equation for ^ having the form 

6^K = K^, [e^^ - en + e^^d^K^, , (48) 

with 5^ K as given by (^) and where we have introduced the recalibrated dimensionless material 
stress-energy contribution 

e'^^ = :^r^^ (49) 

oc 

When e^jy = and e^v = 0, equation (|48| ) can be interpreted, using (^), as the equation of 
evolution for a scalar field ^ non-minimally coupled to the Ricci scalar, i?, and with an effective 
mass given (of the kind discussed by Garriga and Vilenkin |jl6[) by 

m] = ^— [p C^^C^, - (p + \WS^,) . (50) 



4.1 Simulation of Newtonian gravity 

Let us now consider the field in the neighbourhood of a concentrated non-relativistic source, 
that is one whose energy-momentum contribution is given in terms of a preferred timelike unit 
reference vector by the approximate formula 

t'"' -ff"" c^dpu^'u'' . (51) 

We consider a configuration that deviates from an almost uniform and low curvature reference 
configuration by the presence of a matter distribution within a region which is small compared 
to the reference curvature scale, so that we can work at first order both in Cfj^jy and in the 
displacement ^. Hence, the second term of the right hand side of ([48| ) is of order 0{e^) and will 
thus be negligible compared with the first term. On short enough lengthscales, we can keep only 
the gradient terms of highest order, so that we have that d^K and the preceding source 

equation for ^ will reduce to the simple form 

m^^Sp, (52) 

oo 

using the notation 

□ = V'^V;, , K,, = K^.u^'Jf . (53) 

This approximation is the same as the one used when one neglects the cosmological fluids with 
respect to the mass of the Sun on length scales of order of the Solar system. 

The operator □ is interpretable as the wave operator of the local spacetime geometry on the 
brane, and Kqq is the time component of the extrinsic curvature with respect to the preferred 
rest frame, which is interpretable simply as the component of its acceleration orthogonal to 
the worldsheet. In terms of p-dimensional worldsheet coordinates u*, with respect to which the 
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induced metric has components jij and determinant ||7||, (|53| ) can be replaced by the more 
exphcit expressions 

□e = hir'/' (h^^^l''^^) , , ^oo = ■ (54) 



Using (^), (P) and (|T3|), the resulting variation of the first fundamental tensor 7^^ can be 
seen to be 

S^j^, = -2K^,C- (55) 

Since, on the local scale under consideration, the unperturbed intrinsic metric can be taken to 
be of flat Minkowski form, = r]ij say, the corresponding perturbed form of the induced metric 
will be expressible in the form 

7ii = Vij + hij , (56) 
in which, by (^), the induced metric perturbation will be given by 

h,j = -2Kij^. (57) 

For a spherical or hyper spherical static source distribution with integrated mass M, the 
equation (p2[) for ^ has the well known solution 



^00 \ M 



-p-3 ' 



(58) 



for p > 4, while for the case of a circularly symmetric source distribution on a membrane in an 
ordinary 4-dimensional spacetime background the solution is 

^- ^°°-Mlnr, (59) 



where, in both cases, r is the radial distance from the centre. 
The rest frame component 

= hiju'u^ , (60) 

of this perturbation is given in terms of the corresponding rest frame component Kqq of the 
second fundamental form by the dimensionally generalised Newtonian limit formula 



M 

Ko = 2G[p] =^ , (61) 



if one identifies the unrationalised p-dimensional gravitational constant as 



K. 



2 



We now concentrate on the constraints that have to be taken into account in the case p = 4 
which is of greatest interest to us but we will also we then discuss some applications of the p = 3 
case. 
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4.2 Observational constraints when p = A 



It seems from (|60| ) and (q^) that this mechanism can provide a plausible theory for Newtonian 
gravitation but there are some important additional constraints that have to be taken into 
account. In this section we shall start by considering the constraint arising from the requirement 
that as well as reproducing Newtonian effects, the theory should also have the correct relativistic 
behaviour in the weak- field limit, meaning that it should for instance agree with solar system 
measurements of gravity. 

With respect to an unperturbed metric of the standard Minkowski form, 7]^^ = —1, rjQa = 0, 
Vab = ^ah (using Latin space indices, a = 1,2,3) a perturbation hij of the form introduced 
in ( ^6|) will be expressible in terms of the relevant parameters of standard post-Newtonian 
parameterisation [|l^ by the formulae 

h,, ~ 2[/-2/3C/2, 

hoa --{ai - 2a2)waU , 

Kb ~ 2-fU5ab (63) 

so that the trace h = h\ will be given by the formula 

/i = 2(37 -1)[/ + 2/3 (64) 

in which the quantities Wa are the velocity components of the PPN coordinate system relative 
to the mean rest-frame of the universe (which is usually interpreted to be that of the cosmic 
microwave background) and U is the Newtonian potential, whose asymptotic form at large 
distances from the source will of course be given by [/ ~ Gj^jM/f. The standard post-Newtonian 
parameters 7, /3, ai and 0:2 are fixed physical constants. 

Any gravitational theory should predict values for these constants and confront the results 
with what is observed. In particular, Einstein's general relativity gives that 7 = 1, /? = 1 and 
that all other post-Newtonian parameters must be zero. Less restrictively, Jordan - Brans - 
Dicke theory involves a free parameter u in terms of which 7 = (cj + l)/(a; + 2), while giving 
(3 = 1 and zero values for the other post-Newtonian parameters |jl^. The freedom to adjust 
the parameter oj can be used to match the the Jordan-Brans-Dicke theory to a generic 



configuration of the form (57). However this freedom ceases to exist when measurements of 
gravitational light defiection, and the Shapiro time delay effect, are taken into account. These 
effects provide an observational bound on the PPN parameter 7 whose current value is given by 

I7- 1| < 3 X lO"''. (65) 

This quantity measures how much space-curvature is produced by a unit rest mass, whereas 
the quantity (3 measures the degree of "non-linearity" in the superposition law of gravity, an 
effect whose estimation is beyond the scope of the purely linear analysis provided here. By 
substituting ( |57D in (^) and ( |64[ ) it can be seen that what our linear analysis does provide is a 
limit 

K 



< 9 X 10"^ (66) 



that must somewhow be satisfied by the extrinsic curvature as a condition for compatibility with 
the observational limit ( |65[ ) on the assumption that the unperturbed configuration is spatially 
isotropic. 
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Apparent deviations from isotropy might arise from the effects of the parameters ai and 02, 



which (due mainly to Lunar ranging experiments) are subject to constraints [17| given by 



ai < 2 X 10-"^ , 02 < 4 X 10"'^ . (67) 

In order to be compatible with these limits also, in addition to satisfying the limit (|6^), the 
unperturbed extrinsic curvature must satisfy the restriction 

^ « -\a,wa (68) 

with ai subject to ( |67| ) where the velocity vector Wa has magnitude of the order of 200 kms~^ 
while its direction is opposite to that of the cosmic microwave background relative to the Earth. 

Before considering the further restrictions imposed by global cosmological considerations on 
this model for effective gravity in a brane-world scenario, we shall briefly digress to discuss 
analogous local effects in 2-branes. 

4.3 Application to the case p = 3 

As a second application, consider the well known "museum" experiment used to illustrate how 
gravity is a consequence of spacetime geometry and in which a heavy bead of mass M is placed 
on a piece of fabric of tension, 1^ say, in the local terrestrial gravitational acceleration field with 
magnitude g. The weight of the heavy bead produces a dip in the level of the surrounding fabric, 
and thus creates an effective potential well in which another lighter bead may be observed to 
orbit the central point of attraction in much the same way as a planet orbits the Sun, but with 
an acceleration that varies inversely as the distance, not inversely as the square of the distance 
as in the Keplerian planetary case. In this example, at lowest order the component of the force 
density orthogonal to the worldsheet (vertical in this case) is 

7 = X.9, (69) 
and the corresponding displacement will be given by 

e=^Mlnr, (70) 



so that the effective gravitational potential is 

I 
27rT 



y.« = l^Mlnr, (71) 



as long as we are sufficiently far from the centre for the weak-field regime description to be 
applicable. On approaching the centre, on enters a strong-field regime. 

A less trivial but much more speculative application is implicitly suggested by the recent 
proposal [^] that the rotation curve of spiral galaxies may be explained by supposing that these 
galaxies are trapped on lightweight domain wall, of which the thickness is comparable with the 
thickness of the luminous part of the galaxy. The idea of Dvali et al is that the domain wall 
supports a massless scalar particle zero mode confined to (2-|-l)-dimensions with a coupling to 
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the galactic matter so that at distances larger than the typical thickness of the wall the total 
potential for the net force that acts will have an extra logarithmic contribution. The formula 



v'=r^. (72) 
dr 

can be used to obtain an observational estimate of the radial dependence of the effective gravi- 
tational potential V from measurements of the galactic rotation velocity v. In typical cases the 
result can be described in terms of a plateau at distances larger than the radius of the galaxy 
disk where v tends towards a constant limit so that one obtains an empirical formula of the 
form 

V=-^Mo^+vllnr, (73) 

with V ranging from 10~^ to 10^'^ in the relativistic units we are using, that is ranges from 
60 kms-i to 300 kms'^ It was suggested in pOf that these observations be used to tune 



the parameters, such as the energy scale, of their light scalar field, which does not imply any 
cosmological catastrophe. 

What we wish to point out here is that this scalar field mechanism is very similar in its 
effect to the simulated gravity mechanism we have been describing here. According to our 
preceding analysis, if the galaxies were constrained to move on a 2-brane subject to an accel- 
eration corresponding to a local extrinsic curvature component Kgg then, in addition to the 
genuinely gravitational inverse square law contribution, there would be an extra contribution to 
the effective gravitational field so that the effective gravitational potential would be given by 



Kff = -^ + -;^M^^,lnr, (74) 



This would evidently have the same effect as the potential (|7^) if the acceleration and tension 
of the membrane satisfy 

-'oo -"^Gal 

Of course, in order for this variant of the original mechanism to be made astrophysically plausi- 
ble, it would still be necessary to find some natural mechanism for actually trapping the galaxies 
on the 2-brane. 



5 Cosmological matching conditions. 

5.1 Effective Friedmann equation 

It should be remarked that none of the preceding results depend on the specific functional form 
of the background metric coefficient V in (23). The most general solution, as derived in |^] and 
discussed in the preceding work [Q, ^] is 

V-/, 2^ ,2 2Gt,+i].M ^^^^ 

p{p — 1) 2 

The first term is the curvature, k, of the homogeneous isotropic space metric di'^, the second term 
is the effect of the cosmological constant. A, of the bulk and the third term is the contribution 
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of the total asymptotic mass, A4. Substituting this in (|25D, and using the solution (pq) of the 
worldsheet evolution equation (^), we immediately obtain the relevant Friedmann type equation 
for the comoving lengthscale r = a of the brane in the form 

«V 1 (7 , ^bV , 2A 2G^^^M k 
a) W\p aP) p{p-l) aP ' ^ ' 

As observational evidence suggests that the spatial curvature of the universe does not play a 
significant role, we shall assume that A; = from now on. When the gravitational coupling 
is very small, the asymptotic mass would have to be very large to have a significant effect so 
we shall also set = 0. However, in the case where there is no gravity, this term could be 
reinstated with a different physical interpretation. 

In order to match our model to what is observationally known about the large-scale structures 
of our actual universe, an obvious first step is to set the brane worldsheet dimension to p = 4. 
As in the preceding discussion Q of the gravitationally coupled case, the low density epoch in 
which we exist can be dealt with by an expansion in powers of the dimensionless density ratio 

e = ^. (78) 



Equation (^) can be rewritten as 



r2 




(79) 




where the effective 4-dimensional cosmological constant, A^, is given in terms of its 5-dimensional 
analogue. A, and the force, /, by 

- A / T \2 

(80) 

It can immediately be verified that (|8^) agrees with what is obtained from the corresponding 
formula in the preceding analysis when the background gravitational coupling constant is set 
to zero. However, as we have already pointed out, if we similarly set the background gravitational 
coupling to zero in the corresponding formula Q for the effective 4-dimensional gravitational 
coupling constant G^^j we would obtain a negative value that cannot possibly matched to the 
Newton constant Gj^ that is actually observed. This is because the scheme of the previous work 
corresponds to assuming that the constant is sufficiently small for its contribution to be 
entirely neglected in the late time limit considered here, so that the bracketed contribution in 
(|79| ) contains at lowest order just the £ term, which has the desired property of proportionality 
to p, but with the wrong sign. 

The idea of the present approach is to get round this apparent obstacle to taking the zero 
background gravity limit by making an essentially different hypothesis, namely that instead of 
being small, is large enough to give a positive contribution that can overcome the negative 
contribution proportional to —£. The cosmological scenario will thus have an era dominated by 
what we refer to as pseudo-radiation, that will last until at least the present epoch the unless 
it was superseded by a transition (which if it occurred must have been recent) to an epoch of 
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domination by a cosmological constant. Instead of the term "mirage matter" that was used 
by Kehagias and Kiritsis ||5| whose pioneering analysis of gravitationally uncoupled brane-world 
scenarios originally drew attention to the potential cosmological importance of such an effect, we 
prefer to use the term "pseudo-radiation" for the background geometric contribution dominating 
the effective Friedmann equation because, like a true radiation gas, its energy density evolves 
proportionally to a~^, while unlike an ordinary mirage its dynamical effect is very real - in our 
case actually preponderant. 



5.2 The nucleosynthesis condition 



In order for the kind of scenario described in the preceding subsection to be realistic, we need 
at least to require that the four parameters {A,T^, f , E^) are such that at nucleosynthesis, the 
pseudo-radiation term dominates the Friedmann equation (in particular compared with terms 
in a~^) and that its amplitude has the value required to get the correct nuclear abundances, 
i.e. when the temperature, B, is comparable with the nuclear reaction value, ©j, 
expansion rate must be given approximately by 

'd\2 SvrG, 
.a. 



-P. 



1 MeV, the 



(81) 



where is the radiation density as given by the formula 

^2 



TT 

30 



(82) 



in which the value of the number g^, of effective relativistic degrees of freedom is calculated to 
allow for the photon field and three families of degenerate neutrinos which gives = 10.75. 
In order to avoid getting a negative value for the right hand side of (^9|) at some subsequent 
time, the parameters must be such that the pseudo radiation contribution proportional to a~^, 
or equivalently to G^, must be larger than the true radiation contribution given by (^) and also 
than the matter contribution in the decomposition p = p^+ p^ In order to be consistent with 
the standard picture of development of clustering from small perturbations the cosmological 
constant must remain subdominant until at least the present epoch. 

To determine the implication from these requirements, we expand ( [79| ) to get 

7^B , K ^Jj E^'''^ 



2 T 



3 T2 



(83) 



and make the decomposition e = + e^, where = p^jT^ and = pJT^- These terms 
will be given in terms of their present values by = ^J^??^ and = e^i)'^, in terms of the 
reduced temperature, as defined in terms of its present value ©o or the present value oq of the 
cosmological scale factor by 

*=^ = =^. (84) 

Bo a 

is thereby expressible in the form, after setting ao = 1, 



The equation (|8 

2 




^^ + E^(E^- fe 



1} 



(85) 
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where we have introduced the matter-radiation density equahty value 



K^^, (86) 



of the reduced temperature, which is much larger than unity (t?^^ ~ 10^ in standard cosmology) 
and in terms of which the total matter component takes the form 

e = eyH^ + ^J- (87) 

Note that the terms of 0{e'^) in (^), although containing terms of order -d^, are suppressed by 
an additional factor of and thus can be dropped legitimately (see below). 

In order for the term proportional to t?^ to remain positive until at least the present time, 
1? ~ 1, we evidently need that the pressureless matter term be smaller that the pseudo radiation 
term proportional to , giving the constraint 

AE, > Je%^ (88) 

and furthermore, to avoid premature domination by the cosmological constant, we also need 

7{AE^-7eyj>l\A,\T^. (89) 

This implies that the term in -i?^, which has the wrong sign, must always have been subdominant 
until now: it can at most equal the pseudo-radiation component today, and thus was much less 
in the past since it redshifts more slowly. 

Since the true radiation term is always a fixed fraction of the pseudo-radiation term, the fact 
that our scenario is compatible with the nucleosynthesis data leads first to the identification of 
the Newton constant, Gj^ in the formula (^) as 

Since 'd^^ ^ 1, the positivity of Gj^ is evidently ensured by (|88|). It is to be observed that 
this globally determined quantity Gj^ is a priori different from the quantity Gj^j determined 
according to (^) by local physics, so the requirement that they should ultimately agree will 
entail a rather excessive amount of fine tuning. 

In order for the term in (^) proportional to t?^ to be negligible compared to the one pro- 
portional to i?^ at the epoch of nucleosynthesis we evidently need to have 

J{4E, - Je^) > 8E, [E, - /e^X (91) 

where ~ lO^'' is the value of d at nucleosynthesis. In conjunction with (^), this latter 
constraint simply gives the requirement that the four parameters -Eg, /, and A should satisfy 
the conditions expressible in terms of the present cosmological radiation density 

P° = = |^5*eo\ (92) 
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as the inequalities 

» 2^ > ^ . (93) 

J CG 

The only other restrictions on the four listed parameters are the conditions obtained from the 
inequality (^9|) and the identification (^) — in which the second term on the right is negligible 
by (|8^) — so that we are left with 

— ~ - — - > -A 94) 

P oo 

using the standard definition, nip = G^^^'^, of the Planck mass, whose numerical value is given 
by nip ~ 10^2 MeV. No further restriction is needed to ensure that the term proportional to ■i?^^ 
in (|83| ) remains negligible at temperatures up to the nucleosynthesis range. 

It can now be seen to follow directly from ( |9^ ) that the brane mass scale m^, defined by 
= 7^ , must satisfy 

m^^{p%,y^'^^- (95) 

With the standard values t?^^ ~ 10'^, i}^ ~ 10^'^ and Bq ~ 2 x 10^^'^ MeV, which corresponds to 
~ 6 X 10-39 (MeV)^ we get that 

> 10 MeV (96) 

which is comparable with what we found in our previous work Q. The bound (|9^ ) implies that 
<C 10"^'* which justify a posteriori that we can neglect terms of order in (^). 
The first equality in (^) will now give an order of magnitude inequality for the product 
fEp , which will be given in terms of the brane tension by 

f^B in-81/A/r„-lA2 



T2 



IQ-'^'iMeVy, (97) 



and hence by ( |96[ ) will satisfy fE^ > 10~^3(MeV)^'^. In order for the cosmological constant to 
have become important by the present epoch (as suggested by some recent supernova data) we 
would need that the ratio A^/Hq should be comparable with unity where Hq is the present day 
value of the Hubble expansion rate which is given roughly by where Hq ~ 10"^^ MeV. However 
it follows from ( |9^ that this condition would fail by a considerable margin since substitution of 
& gives 

AJH^ < 10'" . (98) 

The limits on Ep obtained from ( |93[ ) by substitution from the first equality in (94) have the 
form 



IO'^M^TFTt) MeV^ > > IQ-^M MeV^ (99) 



10 MeV / " V 10 MeV 



2 



6 Discussion 



It seems from the former paragraph that we can find a set of parameters for our model to be 
compatible with large-scale cosmology (that is, with nucleosynthesis and a cosmological constant 
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dominated universe today), even if very contrived. But it has to be emphasized that the large- 
scale value of the gravitational constant 



3 / jE^ 
87rT3 2 i eO 



7 

4 



(100) 



has no reason a priori to agree with its small scale value (32) 

G 



[4] 



47rT. 



(101) 



We first investigate here, the compatibility between the local and global properties of the 
gravitation and then discuss the viability of such a model as well as its possible extensions. 



6.1 Compatibility between local and global matching 

To compare the two values of the Newton constant, respectively determined from the local 
Newton law and global cosmology, we need to compute the extrinsic curvature in the cosmological 



framework. From (38) we get, at first order in e, that the extrinsic curvature scalar is given by 



i^(global) 



+ 



6 + r 



(102) 



on the assumption that, as described in more detail in the appendix, one assumes a polytropic 
equation of state with index F for which the pressure of the cosmic fluid in ( pT[ ) will be given by 
P = (F — l)yO. Equation (p^) then implies that we shall have 



jv^{global) 

-fi-nn 



4T 



/-3(F + l)U + i?3^4 



(103) 



In the relevant range (|93|) for the admissible parammeter values, the large-scale value of the 
extrinsic curvature components will therefore be given at lowest order in e today (i.e. when 
i9 ~ 1) by 



(global) r-^ _J_ 



f 7^^(global) ^_ / 



(104) 



The extrinsic curvature vector will thus be subject to what is interpretable as a space-time 
isotropy (Lorentz invariance) condition, namely 



T^(global) 



4 



(105) 



which is evidently quite different from the observational constraint (^6|). However the apparent 
contradiction between these two conditions need not be of any consequence if, as we have already 
concluded, the relevant local and global curvature values are independent. 



To reconcile the cosmologically determined value Gj^ given by (100) with the locally measured 
value Gj^j, given by (101), of the gravitational constant we need to have 



7V'(local) 



4T2 eO 



1/2 



(106) 
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with the supplementary requirement 



2^(local)_ 



This local value is obviously different 



from the value (104) found on large scales since we shall have 



J^(global) 



a 



T^(local) 

-"■00 



with EE 12e° ^ 



(107) 



where by (pq), we get that the parameter a is in the range 



2 X 10^3 > > 9 X 10^3 



10 Mev / 



(108) 



from which it follows, a fortiorei, that a < 1. 

Such mismatch between the small-scale and large-scale values of the extrinsic curvature does 
not seem unreasonable in view the consideration that the background configuration of brane 
does not need to be strictly homogeneous and isotropic, as assumed until now, but may have 
different curvature on small and large scales. The cosmological constraints determine the large- 
scale extrinsic curvature, while light deflection observations so far give precise information only 
about local extrinsic curvature. This raises the question of the transition scale between these two 
regimes which should range somehow above the galactic scale. This also leads to a potentially 
lethal test for our model: since the local extrinsic curvature has no reason to be the same in 
two different galaxies, the measurement of the light deflection should lead to different results in 
different galaxies. 



6.2 Viability of the model 

We have shown that as well as automatically providing a local gravitational interaction of the 
familiar Newtonian kind, the simulated gravity model presented here can be contrived in such 
a way as to satisfy the most stringent cosmological constraint - namely the one provided by 
nucleosynthesis. However it is still confronted with other serious problems. 

• Firstly, it suffers from a double fine tuning requirement arising from the fact that, on local 
scales, the gravitational constant depends a priori on time. We must require not only 
that the weak-field limit be instantaneously compatible with the the condition ( |66[ ) for 
consistency with the observed deflection of light, but also that the predicted gravitational 
constant (^) should vary slowly enough to be consistent with the experimental limit [23 



G 



[4] 



< 6 X lO^^Vr"^ 



(109) 



Secondly, the predicted post-Newtonian effects and the current observational limits [17 
imply that the local value extrinsic curvature is subject not just to the limitation its 
mixed components K^^a given by ( |67| ) and (65) but on the further fine tuning condition 
subjecting its time component K^^ to (p6|), which is eqivalent to the condition that the 



(110) 



extrinsic curvature should turn out to be almost exactly of the form 



entailing v' ~ —r'ja in (p^). Note that, at this stage of our investigation, these restrictions 
are not worse than those occuring e.g. for scalar-tensor theories of gravity (see [jl^ for 
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a comparison) before the mechanism of attraction |21| toward General Relativity was 
discovered. However at this stage we are not aware of any natural mechanism that would 
lead to the form ( |110D satisfying the post-Newtonian constraints as a preferred outcome. 

In addition to the kinds of observational restriction we have considered so far, it is to be 
noted that there will be others involving non-linear effects and motion of the source. It is 
not our purpose to investigate all these constraints here, but just to point out that they are 
of potential importance whenever alternative theories of gravity are under consideration, 
and that merely recovering the standard Newton law is in general not enough. None of 
these post-Newtonian constraints have been explicitly taken into account so far in this 
or other works on brane-world scenarios. However, in typical brane-world models, (even 
those of asymmetric type) gravity at Solar system scales is describable in terms of a spin 
2 graviton field, and in such cases it can be expected that the post-Newtonian constraint 
will automatically be satisfied. 

Thirdly, concerning the cosmology of our model, the expansion rate switches from a pseudo- 
radiation to a cosmological constant dominated phase whereas the matter content switches 
from a radiation dominated era to a matter dominated era (at i?^,^) and then to a cosmo- 
logical constant dominated era (at i? 1). This is a non-standard cosmological scenario 
that does not have an ordinary matter-dominated era characterised by a Friedmann type 
evolution equation of the form oc a^. raises the question of how large-scale structure 
will develop in this kind of scenario. In the case of the standard cosmology, it is known 
that density perturbations needed for subsequent galaxy formation cannot grow during 
a radiation dominated era. However it is not yet clear what may happen in a pseudo 
radiation era in which the strength of gravity on shorter scales is affected by the nature of 
fluctuations on longer scales. 

Fourthly, in relation to the preceeding remark, we need a transition scale between the local 
and global behaviour that has to range a somewhere above the galaxy scale. However, this 
scale is not directly determined by the model and might involve a third fine tuning. A 
potentially relevant test that can be thought of is weak gravitational lensing on cosmolog- 



ical scales 1 22 1 which can provide a probe of the Newtonian law of gravitation (and more 
particularly of the Poisson equation) up to some hundreds of Megaparsec when data are 
available. This may be able to provide constraints on this transition scale. 



6.3 Conclusions 

In this article, we have considered a simulated gravity model when the reflection symmetry is 
broken in brane-world scenarios for which genuine (bulk) gravitational coupling is absent, but 
in which the possible effect of a gauge four form coupling is allowed for. A previous investiga- 
tion by Kehagias and Kiritsis |^ of this gravitationally decoupled "probe" limitQ is noteworthy 
for having drawn attention to the dynamical importance of what they referred to as "mirage" 
matter, meaning the contribution we have preferred to describe as pseudo-radiation. Their work 
was restricted to homogeneous isotropic cosmological configurations of any codimension, and 
allowed only for the kinds of matter field most directly motivated by higher dimensional su- 
perstring theory, namely a Born-Infeld type contribution in the brane and a Ramond-Ramond 

The motion of these they sometimes describe as "geodesic motion of a probe D3-brane" , which is potentiaUy 
misleading as the evolution is not generated by geodesies in the usual sense of the term. 
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gauge field outside. In this investigation we restrict attention to the codimension one case but 
progress further in two ways. The first advance is by allowing for more realistic matter fields on 
the brane, and showing that the so-called "mirage" contribution, which we have called pseudo- 
radiation, can be adjusted to give consistency with nucleosynthesis — the most empirically 
sensitive cosmological constraint. The second kind of progress made here is the less mathemat- 
ically trivial step of allowing for local deviations from homogeneity and the demonstration that 
in the weak- field, small-scale limit there will be an effective "simulated" gravity mechanism with 
a Newtonian type inverse square law behaviour. 

The outcome of the first advance is rather satisfactory: having derived the equation of 
evolution of the universe we show that, for a certain range of parameters, it can be made 
compatible with nucleosynthesis subject to the proviso that the cosmological constant is not yet 
dominant today. However the outcome of the second kind of progress raises as many problems as 
it solves: as discussed in the previous section, agreement at the local level can only be obtained by 
adjusting the relevant parameters in an artificially contrived manner. It necessitates at least two 
fine tunings (arising from the time variation of the gravitational constant and of the compatibility 
with the post-Newtonian limit), and requires an undetermined transition lengthscale between 
local and global properties of our 3-brane. Let us emphasize that it is not enough to recover 
Newton gravity on local scales and to produce a viable cosmology; one also needs to take the 
post-Newtonian constraints into account. Further work would be needed to derive the detailed 
evolution of density fiuctuations in our model, but it does not seem that it would lead to an 
acceptable description of structure formation, at least if the standard picture is still valid (which 
is not entirely obvious). 

Different, more sophisticated, scenarios might be conceived, in which one might try to over- 
come the problems that have been pointed out by allowing more general fields to act in the bulk. 
Another possibility is to change the background geometry, which we can do without changing 
the bulk matter content when there is no gravity. By tuning the background, it is very likely 
that we could get a cosmology with the standard Friedmann law. In any case, even if the ef- 
fective gravity mechanism described here is not dominant, it will always be present to some 
extent as a contributing effect in refiection symmetry breaking brane- world scenarios. Whereas 
most previous work has concentrated on strongly self gravitating reflection symmetric case, the 
present article initiates investigation of the opposite extreme limiting case of brane-worlds that 
are very weakly self gravitating but for which gravity is simulated by the effect of refiection 
symmetry violation. 
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A Variational formulation in static background 

The variational principle for a minimally coupled 3-brane in a given background specified by 
a metric Qf^i, and a (Ramond-Ramond type p) gauge 4 form A^^p^j will be given |^] in terms 
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of a Lagrangian density L depending only on internal fields on the brane worldsheet, and of a 
"brawn" coupling (generalised charge) constant e^^j , by an action integral of the form 

(£ + ^^M-PM^,,.)d5, (111) 

where d5 is the element of worldsheet volume measure as given in terms of internal coordinates 
cr*, (i=0,l,2,3) and the determinant I7I of the induced metric 7^^- by 

AS = II7II dcj" dcji dcj2 d(j3 , (112) 

and where £f^'^P'^ is the corresponding antisymmetric unit worldsheet tangent tensor (normalised 
so that £'^^^" S^upa = 4!) as specified in terms of the associated internal measure tensor Ehijk by 

^pupa ^ ghijk ^p^ ^a^ _ (^-^3) 

For a bulk metric of the static form ( |23| ) the space metric ds^ = d£^ of the homogeneous 
isotropic space sections will similarly determine an analogous antisymmetric tensor £i^'^p that is 
tangent to the homogeneous isotropic space 3-sections with normalisation S^'^^Spup = 6. This 
space (as opposed to space time) measure tensor can be used to define a conformally related 
tensor 

ypup = r£f,up, (114) 

with the noteworthy property of having a covariant derivative that is entirely antisymmetric, 
which qualifies it for description as a Killing Yano tensor. This special property can be seen to 
follow just from the "warped product" form (^) of the full spacetime metric ds^ (independently 
of the special homogeneity and isotropy property postulated for d£^) which enables us to evaluate 
the covariant derivative as 

4 

"^pVupa = -ripyupa] ■ (115) 

It follows that the symmetric product tensor 

JCpu = lyppayr, (116) 

with trace = 3r^, will be an ordinary Killing tensor, meaning that its symmetrised derivative 
vanishes, 

V(^/C,p)=0, (117) 
so that for a free particle motion it determines a quadratic constant of motion of the kind 



recently discussed in the general context of "warped product" metrics by Droz- Vincent [24|. 

The utility of the Killing Yano tensor y^up for our present purpose is that it can be used to 
construct an explicit static solution for the "brawn" gauge field Any which may be taken to 



be given in terms of the gradient of the time coordinate t in the expression (23) for the static 
bulk metric by 

Apupa = Ft^^py^pa-^, (118) 

where F is a constant, which will determine the value of the force density / on the right hand 



side of the equation of motion (2C) according to the prescription [Q] 

7 = e[4]F, (119) 
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and in terms of which the corresponding gauge independent field tensor, 



can be seen to be given simply by 



(120) 



where e^yp^r is the background spacetime measure tensor, which will be expressible in terms of 
the coordinates introduced in (|2^ ) by 



'^pvprjT 



20r[^t,,4,,]. (122) 



Let us now specialise, as in § ^, to the case in which the brane is itself homogeneous and 
isotropic, so that the evolution of its worldsheet will be describable simply by specifying r as 
a function of t. In this case the worldsheet tangent tensor £^^^P^ will simply be specifiable in 
terms of the preferred cosmological unit vector u'^ introduced in (p8|) by 



and hence by (|118| ) the coupling term in (111) will be given just by 



^^^^"'^ = ^n'^t,^ . (124) 

In the homogeneous isotropic case (and even more generally for the treatment of adiabatic 
cosmological perturbations) little generality will be lost by supposing that the matter on the 
brane is describable as an irrotational perfect fluid characterised by an equation of state spec- 
ifying the energy density p as a function of the number density, n say of conserved particles, 
which (in the subnuclear energy regime we are concerned with) can taken to be baryons. In such 
a case the 4-dimensional fluid dynamics will have a particularly simple Lagrangian formulation 
in which the only independent variable is a scalar potential, f-p say, with a timelike gradient ^p^p 
whose magnitude /i, as given by 

-p" = -f , (125) 

is to be identified with the relevant chemical potential (i.e. effective mass per baryon) as defined 

by 

= ^ (126) 

while the Lagrangian itself can be taken to be the corresponding pressure function P, as obtained 
in terms of Ji by application of the standard prescription 

P = njL-p. (127) 

To adapt this simple perfect fluid model to the context of a brane-world in a 5-dimensional 
background, all one has to do is to identify the Lagrangian density L in ( |111| ) with the negative 
of the total brane tension T as given by (pi]), i.e. we take 

C = -T, (128) 
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with 

r = T^-P. (129) 

In order to be explicit, for typical cosmological purposes a good description will be provided 
by an equation of state of polytropic form given in terms of a fixed mass (per baryon) m, a fixed 
polytropic index T, and a fixed proportionality constant k say by 

P = + P. , (130) 

with 

p^=mn, p^ = Kff, (131) 

which gives 

fl = m + kTu^-'^ , P = [T - l){p - mn) , (132) 
so that the required pressure function is finally obtained in the form 

, /7i - m\r/(r-i) 
P = ^{r-l){^) . (133) 

The simplest relevant possibility is that of a pure radiation gas, as characterised by m = and 
r = 4/3 for which the pressure function contribution to the Lagrangian density will just be 
proportional to {-^^^ f^ifj)'^ . 

Whatever the form of the equation of state, the four dimensional surface integral in ( |111| ) 
will be proportional in the homogeneous isotropic case to a one dimensional integral of the form 

I = J r^-Tdr + e^^^F^-dt). (134) 



Using the explicit form of the metric (|2^ ) to evaluate the derivative of the proper time r with 
respect to the background metric time t, this will be expressible in the standard form 



J Ldt, (135) 



for an ordinary Lagrangian (not Lagrangian density) function given by 

L^-T[v-V-^Qy''r' + e,„F'^. (136) 

in which the equation of state specifies P, and hence also T, as a function of Jl and hence of 
d^p/dt and dr/dt via the relation 

In terms of the momenta defined by the variation 

the associated Hamiltonian function will be given by an expression of the standard form 

dip dr 
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It is to be remarked that this ordinary kind of Hamiltonian, as defined with respect to the 
background time coordinate t, is to be distinguished from the alternative kind of Hamiltonian 
defined with respect to the proper time r (as used in the Wheeler De Witt type formulation ||^]), 
whose use requires the imposition of an extra constraint (whereby the proper time Hamiltonian 
is initially set to zero in the standard version of the Wheeler De Witt formulation). In the 
present case the relevant momenta are found to be given by the formulae 

P^ = nr\ p^=(v-V-'{^)y^\~'/'Ur^^^, (140) 
in which U is the energy density function given according to (^l|) by 

^ = X.- (141) 
We thus obtain the Hamiltonian in the form 



/ i/drx2\-V2 ^4 

H=[V-V-'{-)j VUr' + e,,,F-. (142) 

The time independence of the background ensures that the motion will be such that this 
Hamiltonian has a constant value 

H = E^, (143) 

that is identifiable as the integration constant introduced in (|3^), which is thus interpretable 
as a cosmological energy density. A more trivial observation is that since the Lagrangian does 
not depend on the undifferentiated scalar (p but only on its gradient (i.e. since ( |138| ) contains 
no dLjdip term) the corresponding momentum "p^^ will also be conserved, a conclusion that can 
be seen from ( |140| ) to be interpretable as equivalent to conservation of baryons, whose number 
density varies inversely as the cosmological volume, i.e. n oc . 

As a final remark concerning energy conservation, it is to be observed that the stationary 



character of the background specified by (^3|) and ( |118|) with respect to the action generated 
by the Killing vector given by kf^dx^ = —Vdt ensures that even when inhomogeneous 
perturbations are taken into account there will still be an energy density current, E say, given 
by 

W = k^ir^^ + ^^^^^^^M-p-) , (144) 

that will satisfy a surface flux conservation law of the simple form 

V^W = {). (145) 
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